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We employ the dynamical mean-field approximation to study the magnetic properties of a model 
relevant for the dilute magnetic semiconductors. Our model includes the spin-orbit coupling on the 
hole bands, the exchange interaction, and the attractive Coulomb potential between the negatively 
charged magnetic ions and the itinerant holes. The inclusion of the Coulomb potential significantly 
renormalizes the exchange coupling and enhances the ferromagnetic transition temperature for a 
wide range of couplings. We also explore the effect of the spin-orbit interaction by using two 
different values of the ratio of the effective masses of the heavy and light holes. We show that in the 
regime of small J^-V the spin-orbit interaction enhances Tc, while for large enough values of Jc-V 
magnetic frustration reduces Tc to values comparable to the previously calculated strong coupling 
limit. 

PACS numbers: 75.50Pp, 75.30.Et, 71.10.Hf, 71.27.+a 



I. INTRODUCTION 

Although the notion of using magnetic semiconduc- 
tors in spintronic devices dates back to the 1960's^, the 
discovery of high temperature ferromagnetism in dilute 
magnetic semiconductors (DMS^P'^ initiated an active 
search for the optimal compound with a magnetic tran- 
sition above room temperature. Since these materials 
are good sources of polarized charge carriers, they may 
form the basis of future spintronic devices,HI51 which uti- 
lize the spin of the carriers as well as their charge to 
simultaneously store and process data. Perhaps one of 
the most promising DMS is GaAs doped with Manganese 
due to its rather high ferromagnetic transition temper- 
ature (Tc > 150 K for bulk samples and ~ 250 K for 
i5-doped heterostructure^^l^ and its wide use in today's 
electronic devices. 

In Gai_a;Mna;As, the Mn+-^ ion primarily replace Ga+'^ 
playing the role of acceptor by introducing an itinerant 
hole to the p-like valence band. The strong spin-orbit 
interaction in the valence band couples the angular mo- 
mentum to the spin of the itinerant hole resulting in total 
spin J— I + s—3/2 for the two upper valence bands and 
J=l — s=l/2 for the split-off band. Each manganese also 
introduces a localized spin (5=5/2) due to its half- filled 
d orbital. In addition, since the Mn"*"^ ion is negatively 
charge with respect to the Ga+^ ionic background there 
is an effective attractive interaction between the Mn ion 
and the charge carriers. 

In previous studie^^ some of us have explored the 
effect of the strong spin-orbit coupling on the ferromag- 
netic transition temperature Tc, the carrier polarization 
as well as the density of states and spectral functions us- 
ing the Dynamical Mean-Field Approximation (DMFA). 
In these studies we used the k-p Hamiltonian to model the 
dispersion of the parent material (GaAs) . While k ■ p is 
a good approximation around the center of the Brillouin 



zone (r point), it is a poor one away from it. In this work 
we improve our model by incorporating a more realistic 
tight binding dispersion for the valence bands as well as 
an attractive on-site potential between the Mn ions and 
the itinerant holes. Moreover, we study the effect of the 
spin-orbit interaction of the holes on the magnetic be- 
havior of the DMS. We find that for intermediate values 
of the exchange coupling both the on-site potential and 
the spin-orbit enhances the critical temperature, while 
in the strong coupling regime the spin-orbit interaction 
significantly suppresses TjlSl. 



The effect of the attractive Coulomb potential has been 
discussed for models with only one valence band, which 
ignore the spin-orbit interaction,!^^!!!! and multi-band 
tight-binding models, which include spin-orbit coupling, 
but with a limited sampling of disorder configurations.'^^ 
Here we include on an equal footing the effect of the at- 
tractive Coulomb potential using a simple Hartree term, 
the exchange between magnetic ions and itinerant holes, 
the spin-orbit coupling, and the disorder within the co- 
herent potential approximation (CPA).'iSHIll Wg investi- 
gate the ferromagnetic transition temperature, the aver- 
age magnetization of the Mn ions, the polarization of the 
holes, and the quasiparticle density of states as function 
of the Coulomb and exchange couplings. First, we use a 
single band model where spin-orbit interaction is ignored 
and carriers have angular momenta J — 1/2. Next, we 
introduce the spin-orbit coupling in a two-band model 
with J = 3/2. By changing the ratio of the masses of the 
light and heavy bands (m;/m/j) we explore the effect of 
spin-orbit coupling. This is the minimal model that qual- 
itatively captures the physics of DMS, however, a more 
realistic approach should incorporate the conduction and 
split-off bands and this will be discuss in future studies. 
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FIG. 1. (color online). Spin-dependent density of states for 
temperature T=0.01t, exchange coupling Jc=2 1 and Coulomb 
potential V = lt,2t and 4t. T=0.01t is below the ferromag- 
netic transition for all values of V. 
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FIG. 2. (color online). Spin-dependent density of states for 
temperature T=0.04t, exchange coupling Jc=5t and various 
values of the Coulomb potential coupling V. 



II. MODEL 

We employ the simplified Hamiltonian proposed by 
Zarand and Jank^^^l -with an additional Coulomb poten- 
tial term: 

H = Ho + JcY.S{R,)■JiR,) + vY,<R^), (1) 

i i 

where Hq includes both electronic dispersion and spin- 
orbit coupling of the holes in the parent compound, Jc is 
the exchange coupling, V the Coulomb strength, S(i?i), 
3{Ri) and n{Ri) are, respectively, the spin of the local- 
ized moment, the total angular momentum density and 
the density of the carriers at random site i. Short range 
direct or superexchange between Mn ions is ignored since 
we are in the dilute limit and we are not including clus- 
tering effects. 

As discussed previously,'^ within the DMFA the 
coarse-grained Green function matrix is: 

GiioJn) = ^ Y.^^LoJ - Hoik) t{iuj,,)]-\ (2) 

k 

where N is the number of k points in the first Brillouin 
zone, /i the chemical potential, and Ho{k) and I](ia;„), 
are matrices representing the band structure of the par- 
ent material and the selfenergy, respectively. The mean 
field function t/o(*w„) = [G~^(iaj„) -I- Yi{iujn)]~'^ is re- 
quired to solve the DMFA impurity problem. At a non- 
magnetic site, the Green function is simply the mean field 
function Gnonii^n) — GoiioJn)- The Green function at 
a magnetic site is Gs(ia;„) = [QQ^{iujn) + Jc^ ■ J + V]^^ 
for a given local spin configuration. 

Next we average Gmag over different spin orientation 
of the local moment. The relatively large magnitude of 
the Mn moment justifies a classical treatment of its spin. 



To get the averag e over the angular distribution we use 
the effective actiorPl^U 

5e//(S) = - ^ \ogdei[go\iujn) + JcS • J + F]e*"""^. 

n 

(3) 

The average over spin configuration is 

(G,„g(iw„)) = ^ J dflsGsiiuJn)exp[-Scs{S)], (4) 

where Z is the partition function, Z = 
J (irisexp(— Scff (S)). Finally the disorder is treated 
in a fashion similar to the coherent phase approximation 
(CPAJP^ESl and the averaged Green function reads 

Gavg{i^n) = {Gmg)x + ^o(«Wn)(l " x) whcrC X IS the 

doping. 

We obtain the hole density of states from the coarse- 
grained Green function in real frequency domain: 

Gm^^J2^^i-Ho{k)-m)r' (5) 

k 

where fl = uj + iO+. The total density of states (DOS) is 
DOS{n) = --ImTrGin), (6) 

TT 

where Tr is the trace. Each diagonal element of the 
Green function ( ImG{^)) corresponds to the density 

TT 

of states for a specific Jz component. 



III. RESULTS 

Since Gai_a;Mna;As is grown using out of equilibrium 
techniques a noticeable fraction of manganese lies not 
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FIG. 3. (color online). Ferromagnetic transition temperature 
versus magnetic exchange coupling for various values of the 
Coulomb potential. Inset: polarization of the holes as func- 
tion of T/Tc for a wide range of values of Jc and V . Notice 
that all the polarization data collapse on a single curve. 



on the Ga site (substitutional) but on the As site (anti- 
site) or somewhere in the middle of the crystal structure 
(interstitial)^^. The real nature of inte rstitia l defects is 
still controversial and yet to be resolved j^^ISH but the one 
consensus is that in most samples there is strong com- 
pensation of the holes introduced by substitutional Mn. 
The density of carriers can also be controlled with elec- 
tric fields.'^ We take these considerations into account 
by simply setting the filling of the holes to half of the 
nominal dopinglSI. We focus on the doping x—5% and 
hole filling rih = x/2. 
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FIG. 4. (color online). Density of states for V—O. and 
Jc =0.5 ti, T=0.005ti; and Jc =0.75 at a temperature 
T=l/130t;. The impurity band is well formed with Jc=0.5ti 
and increasing of the coupling shifts the impurity band to 
higher energies. 



We start by discussing a simplified one-band model 
where we ignore the spin-orbit interaction. Our carrier 
dispersion is e^. = — 2i(cos(kx) -l-cos(ky) -|-cos(kz)), where 
t is the spin independent hopping integral. Fig. [T]and[2] 
display the spin-dependent density of states (DOS) close 
to the edge of the valence band for coupling constant 
Jc—2 1 and 5 1, respectively. Note that inclusion of the 
spin-independent attractive potential results in shifting 
the energy of the holes (electrons) to lower (higher) en- 
ergies for both spin species. This is in agreement with 
previous studieJ^^Hl pjg jY] illustrates the strong influ- 
ence of the Hartee term on the states close to the va- 
lence band edge for moderate exchange coupling. It is 
clear that increasing the Coulomb potential accelerates 
the formation of the impurity band and its splitting from 
the valence band. Fig.[2]shows that for couplings as large 
as Jc=5 1 the impurity band is well formed even for rel- 
atively small Coulomb potentials {V=l t) and the mere 
effect of the Coulomb term is to shift the impurity band. 
Notice also that the predicted shift of the impurity band 
is too large. We believe that this is a consequence of ex- 
cluding the conduction band from our model, since band 
repulsion with the conduction band pushes the impurity 
band to lower energies. 

The main panel in Fig. [3] shows the dependency of Tc 
on the exchange coupling for different Coulomb poten- 
tials within this simplified one-band model. Comparing 
this figure with Fig. [T] and [2] it is clear that Tc increases 
as impurity band forms and separates from the edge of 
the valence band. For each value of V we can identify two 
values of for which the slope of the vs. Jc curve 
changes. For Jc < Jmim Tc increases very slowly, for 
Jmin < Jc < Jsat the impurity band begins to develop 
and Tc increases with the largest slope, for Jc > Jsat the 
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FIG. 5. (color online). Density of states for V=ti, Jc—0.5ti 
T=0.005t;; and Jc =0.75 i; at a temperature T=l/130ti. The 
attractive potential enhances the formation of the impurity 
band as compared with Fig. |4] The chemical potential lies 
in the middle of the first impurity band, as it is displayed by 
the vertical black line. 



4 







1 

- 1 

1 

1 

1 




r J^=U. /jtj, v=u, i=yin) 

a=0.14 |_ _ 50t^^ v=t|, T=0.005tj " 

_ /=0.75t„ V=0, T=t,/300 
a-1 00 •! c 1 1 
" I. . . J =1.50t„ V=t„T=0.01t, . 

c r r 1 


1 

- 1 








\ 1 








■■■•i 








1 








c 


) 


1 


2 3 4 ! 



ro/tj 



FIG. 6. (color online). Density of states for two values of the 
exchange coupling and Coulomb term: Jc = 0.75 1;, = 
and Jc = 1.5 ti, V = ti and temperatures well below the 
ferromagnetic transition, and for a=0.14 and 1.0. 



impurity band is completely split from the valence band 
and the rate of increase in reduces dramatically. In 
brief, the appearance of the impurity band corresponds 
to the large change in the curvature of Tc vs. Jc- Af- 
ter the impurity band is well formed increasing Jc or V 
does not change Tc significantly. In fact, for > 4t 
we can anticipate the saturation of the critical temper- 
ature. This is an artifact of the DMFA and is due to 
the absence of non-local correlations. Inclusion of those 
correlations leads to magnet ic fru stration of the system, 
which in turn suppresses T^IMIll We will come back to 
this point in more detail later when we discuss the two- 
band model. 

Therefore by increasing the attractive Coulomb po- 
tential Tc is significantly enhanced for values of the ex- 
change in a given interval, JminiV) < Jc < Jsat{V), 
where JminiV) and JsatiV) are function of V. This is 
due mostly to the fact that a positive V promote the ap- 
pearance of localized states at the magnetic sites which 
mediate the magnetic order. However, the physics of the 
ferromagnetic state is not modified by V, since the only 
relevant energy scale is given by Tc, as one expects from a 
mean field theory. This is illustrated in the inset of Fig. [3] 
that displays the polarization of the holes as function of 
T/Tc for a wide range of values of Jc and V, showing 
that all the polarization data collapse on a single curve. 
Thus, the effect of V is just to change the nominal value 
of Jc to a larger Jc-^^ ■ 

Now, we introduce a more realistic approach using 
a two-band model. The spin-orbit interaction and the 
crystal fields lift the degeneracy of the p-like valence 
bands into heavy, light and split-off bands. In our model 
we ignore the effect of the split-off band and focus on 
the heavy and light bands which are degenerate at the 
center of the Brillouin zone.^^ Hq is approximated by 
Ho{k) — W {k)e{k)R{k) , where e{k) is a diagonal matrix 



with entries e{k)n,a = — 2t„(cos(kx) -I- cos(ky) -I- cos(kz)), 
with n = l,h the heavy/light band index and R{k), the 
k-p spin 3/2 rotation matrices.l^In GaAs the mass ratio of 
light and heavy holes at the T point is a=mi/mh= 0.l42Sl. 
We compare the results of our simulation for a=0.14 and 
a=l, keeping the bandwidth of the light hole band fixed. 
Furthermore we scale every parameter according to the 
light holes hopping energy (i/), which set the bandwith 
of the hole band. 

Fig. |4] displays the hole density of states close to 
the edge of the valence band for Jc= 0.5 and 0.75 
in absence of the Coulomb potential, and for temper- 
atures well below the ferromagnetic transition tempera- 
ture. One can anticipate that the formation and splitting 
of the impurity band happens for smaller values of Jc/ti 
than in the one-band model. We can explain this by not- 
ing that the total angular momentum of the holes can 
be as large as J=3/2 for heavy holes, leading to a larger 
contribution to the total energy from the the second term 
in Eq. ([!]). Moreover, for a small filling there are more 
available states close to the center of the Brillouin zone in 
the two-band model than in the one-band model. Larger 
number of spin states available to align along the direc- 
tion of the local moment increases the average exchange 
energy and favor ferromagnetism. 

Fig. [5] displays the density of states for the same 
exchange couplings and temperatures, Jc— 0.5 U, 
T=0.005ti, and Jc =0.75 1;, T=l/130t,, but with a finite 
Coulomb potential V— ti . For these values of the param- 
eters a second impurity band appears in the semiconduct- 
ing gap. The appearance of two impurity bands is con- 
sistent with the fact that the model includes two bands 
with Jz — ±3/2, ±1/2. Notice that the second impurity 
band is more populated with light holes ( Jz = 1/2) while 
the first impurity band, with higher energy, is mostly 
made of heavy holes (J^ = 3/2). Since we keep the filling 
of the holes fixed {nh=x/2) the chemical potential sits in 
the middle of the first impurity band, as shown in Fig. [5] 
Thus, as we discussed previoulsy, the shift of the impu- 
rity band will not have noticeable effects on the magnetic 
properties of the DMS. 

To investigate the effect of the spin-orbit interaction 
we introduce a simple toy model which has all the fea- 
tures of our two-band model except that the heavy and 
light bands are degenerate over the whole Brillouin zone. 
Therefore, heavy and light bands have the same disper- 
sion but different total angular momenta jz=±3/2 and 
±1/2, respectively. The different band masses introduce 
magnetic frustratiorP^^ and by setting a^l.O (mh—mi) 
in our model, this magnetic frustration is removed. Since 
ti is fixed, changes in a alters the dispersion of the heavy 
hole band while keeping the light band fixed. 

Fig. [6] displays the total DOS for two values of the 
exchange coupling and Coulomb potential: Jc = 0.75 ti, 
V — and Jc = 1.5 1;, V = ti, and for a=0.14 and 
1.0. Note that for a=0.14 the impurity band is formed 
at lower couplings. Thus, the spin-orbit interaction en- 
hances the formation of the impurity band. We can ex- 
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FIG. 7. (color online). Ferromagnetic transition temperature, 
Tc, vs. exchange coupling Jc, both in units of ti, for different 
values of the Coulomb potential and for a=0.14 and 1.0. 



energy and can be polarized more easily and become 
bonded to the localized moments forming the impurity 
band. For larger values of Jc and V , Jc > 0.81t; for 

V ^ 0, Jc > O.GOti for V = Iti OT Jc > 0.29i; for 

V = 3 1; , the critical temperature for the model with 
a=1.0 surpasses the one for o; — 0.14 in agreement with 
previous findings in the strong coupling regimd^li^] This 
also can be related with the DOS in Fig. |6j where the 
bandwith of the impurity band for a=1.0 is larger than 
for a=0.14. A larger bandwidth corresponds to weaker 
localization of the holes and higher mobility. Therefore, 
they will better mediate the ferromagnetic interaction 
between the magnetic ions and we expect to see higher 
Tc when a = 1.0. For the largest value of Jc and V we 
study Tc{a = 0.14)/Tc(a = 1.) = 0.35 to compare with 
0.48 obtained in the strong coupling limiui^i. 



IV. CONCLUSIONS 



plain this by noting that changing a from 1.0 to 0.14 
decreases the kinetic energy of the heavy holes (with 
becoming more susceptible to align their spin 
parallel to the local moment promoting the formation of 
the impurity band. Fig. |4] and [5] show explictly that the 
heavy holes are the majority of the carriers in the im- 
purity band. On the other hand the bandwidth of the 
impurity band is larger when a=1.0, pointing to less lo- 
calized holes, which better mediate the exchange inter- 
action between magnetic ions. 

Finally we look at the dependence of the critical 
temperature on the parameters of the model: Jc, V and 
a. The results for different values of Jc-V for q;=1.0 
and 0.14 are shown in Fig. [7) Similarly to Fig. |3] we can 
identify for both values of a a range of parameters Jc, 
V where Tc increases strongly. This corresponds to the 
formation and splitting of the impurity band from the 
valence band. For small values of Jc and V, Tc is higher 
for a = 0.14 but as we increase Jc-V the ferromagnetic 
transition temperature for a = 1.0 becomes larger. 
Eventually Tc saturates due to the lack of non-local 
correlations within the DMFA. We can understand the 
higher Tc for a = 0.14 and small Jc,V by looking at 
Fig. |6] For a = 0.14 the impurity band appears at 
smaller values of Jc and V than for a = 1.0. This is due 
to the fact that the heavy holes have a smaller kinetic 



In conclusion, we have calculated densities of states, 
polarizations and ferromagnetic transition temperatures 
for a one-band and two-band models appropriate for 
Gai_2;Mn2;As. We have investigated the effect of adding 
a local Coulomb attractive potential V between the mag- 
netic ions and the charge carriers. The inclusion of a 
Coulomb term leads to the formation of the impurity 
band for smaller magnetic couplings (Jc), in agreement 
with previous studies^ ^ and it significantly enhances 
Tc for a wide range of Jc, without affecting the intrinsic 
physics of the ferromagnetic transition. We also explore 
the effect of the spin-orbit interaction by using a two- 
band model and two different values of the ratio of the 
effective masses of the heavy and light holes. We show 
that in the regime of small Jc-V the spin-orbit interac- 
tion enhances Tc, while for large enough values of Jc-V 
the magnetic frustration induced by the spin-orbit cou- 
pling reduces Tc to values comparable to the previously 
calculated strong coupling limit. 
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